Electron capture rates on nuclei and implications for stellar core collapse 
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Supernova simulations to date have assumed that during core collapse electron captures occur 
dominantly on free protons, while captures on heavy nuclei are Pauli-blocked and are ignored. 
We have calculated rates for electron capture on nuclei with mass numbers A — 65-112 for the 
temperatures and densities appropriate for core collapse. We find that these rates are large enough 
so that, in contrast to previous assumptions, electron capture on nuclei dominates over capture on 
free protons. This leads to significant changes in core collapse simulations. 

PACS numbers: 26.50.+X, 97.60.Bw, 23.40.-s 



At the end of their lives, stars with masses exceed- 
ing roughly 10 M reach a moment in their evolution 
when their iron core provides no further source of nu- 
clear energy generation. At this time, they collapse and, 
if not too massive, bounce and explode in spectacular 
events known as type II or Ib/c supernovae. As the den- 
sity, p, of the star's center increases, electrons become 
more degenerate and their chemical potential /i e grows 
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). For sufficiently high values of the chemical 



potential electrons are captured by nuclei producing neu- 
trinos, which for densities < 10 11 g cm -3 , freely escape 
from the star, removing energy and entropy from the 
core. Thus the entropy stays low during collapse ensur- 
ing that nuclei dominate in the composition over free pro- 
tons and neutrons. During the presupernova stage, i.e. 
for core densities < 10 10 g cm -3 and proton-to-nucleon 
ratios Y e > 0.42, nuclei with A = 55-65 dominate. The 
relevant rates for weak-interaction processes (including 
(3^ decay and electron and positron capture) were first 
estimated by Fuller, Fowler and Newman yij (for nuclei 
with A < 60), considering that at such conditions al- 
lowed (Fermi and Gamow- Teller) transitions dominate. 
The rates have been recently improved based on modern 
data and state-of-the-art many-body models Q, consid- 
ering nuclei with A — 45-65. (This rate set will be 
denoted LMP in the following.) Presupernova models 
utilizing these improved weak rates are presented in 
In collapse simulations, i.e. densities > 10 10 g cm -3 , a 
much simpler description of electron capture on nuclei 
is used. Here the rates are estimated in the spirit of 
the independent particle model (IPM), assuming pure 
Gamow- Teller (GT) transitions and considering only sin- 
gle particle states forproton and neutron numbers be- 
tween Z,N = 20-40 |J. In particular this model results 



in vanishing electron capture rates on nuclei with neu- 
tron numbers larger than N = 40, motivated by the ob- 
servation |5j that, within the IPM, GT transitions are 
Pauli-blocked for nuclei with N > 40 and Z < 40. 

During core collapse, temperatures and densities are 
high enough to ensure that nuclear statistical equilib- 
rium (NSE) is achieved. This means that for sufficiently 
low entropies, the matter composition is dominated by 
the nuclei with the highest binding energy for a given Y e . 
Electron capture reduces Y e , driving the nuclear compo- 
sition to more neutron rich and heavier nuclei, including 
those with N > 40, which dominate the matter compo- 
sition for densities larger than a few 10 10 g cm -3 . As 
a consequence of the model applied in previous collapse 
simulations, electron capture on nuclei ceases at these 
densities and the capture is entirely due to free protons. 
We will show now that the employed model for electron 
capture on nuclei is incorrect, as the Pauli-blocking of 
the GT transitions is overcome by correlations |6j and 
temperature effects 0, Q ■ 

The residual nuclear interaction, beyond the IPM, 
mixes the pf shell with the levels of the sdg shell, in 
particular with the lowest orbital, gg/2- This makes 
the closed gg/2 orbit a magic number in stable nuclei 
(N = 50) and introduces, for example, a very strong 
deformation in the N = Z — 40 nucleus 80 Zr. More- 
over, the description of the B(E2,0 + — > 2+) transition in 
68 Ni requires configurations where more than one neu- 
tron is promoted from the pf shell into the gg m orbit || , 
unblocking the GT transition even in this proton-magic 
N — 40 nucleus. Such a non- vanishing GT strength has 
already been observed for 72 Ge (N = 40) Q and 76 Se 
(N = 42) [13. In addition, during core collapse electron 
capture on the nuclei of interest occurs at temperatures 
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T > 0.8 MeV, which, in the Fermi gas model, corresponds 
to a nuclear excitation energy U ~ AT 2 /8 > 5 MeV; this 
energy is noticeably larger than the splitting of the pf 
and sdg orbitals (E g9/2 ~E pi/2 j 5/2 w 3 MeV). Hence, the 
configuration mixing of sdg and pf orbitals will be rather 
strong in those excited nuclear states of relevance for 
stellar electron capture. Furthermore, the nuclear state 
density at E ~ 5 MeV is already larger than 100/MeV, 
making a state-by-state calculation of the rates impossi- 
ble, but also emphasizing the need for a nuclear model 
which describes the correlation energy scale at the rele- 
vant temperatures appropriately. This model is the Shell 
Model Monte Carlo (SMMC) approach [HI which allows 
the calculations of nuclear properties at finite tempera- 
ture in unprecedentedly large model spaces. To calculate 
electron capture rates for nuclei A = 65-112 we have first 
performed SMMC calculations in the full pf-sdg shell, us- 
ing a residual pairing+quadrupole interaction, which, in 
this model space, reproduces well the collectivity around 
the ./V = Z = 40 region and the observed low-lying spec- 
tra in nuclei like 64 Ni and 64 Ge. From the SMMC cal- 
culations we determined the temperature-dependent oc- 
cupation numbers of the various single-particle orbitals, 
which then became the input in RPA calculations of the 
capture rate, where we considered allowed and forbid- 
den transitions up to multipoles J — 4, including the 
momentum dependence of the operators. This model is 
described in more details in where, however, a smaller 
model space has been used. 

To validate our method at the early collapse condi- 
tions, we have performed diagonalization shell model 
studies for 64,66 Ni, considering the complete (pf) shell 
for 64 Ni, and adopting, for 66 Ni, the (pf) shell for pro- 
tons and the (pf 5/2^9/2) shell for neutrons. We find 
agreement to better than a factor 2 between the present 
(SMMC+RPA) rates and the diagonalization shell model 
rates at stellar conditions (T < 0.8 MeV), for which the 
latter can still be evaluated. 

For all studied nuclei we find neutron holes in the (pf) 
shell and, for Z > 30, non-negligible proton occupation 
numbers for the sdg orbitals. This unblocks the GT tran- 
sitions and leads to sizable electron capture rates. Fig.Q] 
compares the electron capture rates for free protons and 
selected nuclei along a core collapse trajectory, as taken 
from ^3] ■ Depending on their proton-to-nucleon ratio Y e 
and their Q- values, these nuclei are abundant at different 
stages of the collapse. For all nuclei, the rates are domi- 
nated by GT transitions at low densities, while forbidden 
transitions contribute sizably for > 10 11 g cm~ 3 . The 
electron chemical potential fi e and the reaction Q- value 
are the two important energy scales of the capture pro- 
cess. At a given density, i.e. constant /i e , the rate is gen- 
erally larger for nuclei with smaller Q- values. The rate 
is sensitive to the GT strength distribution, if fi e < Q. 
However, /z e increases much faster with density than the 
Q-values of the abundant nuclei. As a consequence the 




1-4 r I I I 1 1 I I I 1 1 I I I 1 1 I I I 1 1 I I I 1 1 I I I 1 1 I , , , 1 

5 10 15 20 25 30 35 4( 
\L e (MeV) 



FIG. 1: Comparison of the electron capture rates on free pro- 
tons and selected nuclei as function of the electron chemical 
potential along a stellar collapse trajectory taken from |T3l . 
Neutrino blocking of the phase space is not included in the 
calculation of the rates. 
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FIG. 2: Electron capture rates on nuclei as function of Q- 
value for 3 different stellar conditions. Temperature and elec- 
tron chemical potential are measured in MeV. The solid lines 
represent the approximate Q-dependence of the rates as de- 
fined in Eq. Neutrino blocking of the phase space is not 
included in the calculation of the rates, pw measures the 
density in units of 10 11 g cm -3 . 



capture rates on nuclei become quite similar at larger 
densities, say > 10 11 g cm" 3 , depending now basically 
only on the total GT strength, but not its detailed dis- 
tribution. This is demonstrated in Fig.[5]which shows our 
calculated capture rates as function of Q-value at 3 dif- 
ferent stellar conditions. The Q- value dependence of the 
capture rate for a transition from a parent state at excita- 
tion energy Ei to a daughter state at Ef (AE = Ef —Ei) 
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is well approximated by |13J 

(1) 

where x = (Q - AE)/T, r\ = (/j, e + Q- AE)/T, K = 
6146 s and B represents a typical (Gamow- Teller plus 
forbidden) matrix element. The quantities Fk are the 
relativistic Fermi integrals of order k. 

At (pY e = 7 x 10 9 g cm" 3 , T = 0.93 MeV), we observe 
some scatter of the calculated rates around the mean (in- 
dependence indicating that several parent and daughter 
states with different transition strengths contribute. For 
nuclei with large \Q\ -values the simple parametrization 
breaks down. However, at pi%Y e > 4, the electron chem- 
ical potential has increased sufficiently that the rates be- 
come virtually independent of the strength distribution 
and are well represented by the average Q-value depen- 
dence JTJ with B = 4.6 and AE = 2.5 MeV. Such a 
parametrization could then be adopted in core collapse 
simulations for even higher densities, when nuclei heav- 
ier than the ones included in the present study start to 
dominate in the composition. 
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FIG. 3: The reaction rates for electron capture on protons 
(thin line) and nuclei (thick line) are compared as a function 
of electron chemical potential along a stellar collapse trajec- 
tory taken from _1.2|. The insert shows the related average 
energy of the neutrinos emitted by capture on nuclei and pro- 
tons. The results for nuclei are averaged over the full nuclear 
composition (see text). Neutrino blocking of the phase space 
is not included in the calculation of the rates. 

Simulations of core collapse require reaction rates for 
electron capture on protons, R p — Y p X p , and nuclei 
Rh = YiXi (where the sum runs over all the nu- 
clei present and Yi denotes the number abundance of a 
given species), over wide ranges in density and temper- 
ature. While R p is readily derived from the calcu- 
lation of Rh requires knowledge of the nuclear composi- 
tion, in addition to the electron capture rates described 
earlier. The information about the nuclear composition 



provided by the commonly used Lattimer-Swesty equa- 
tion of state 01 ' the total abundance of heavy nuclei and 
the average Z and A, is not sufficiently detailed to make 
adequate use of these new reaction rates. Therefore a 
Saha-like NSE is used to calculate the needed abundances 
of individual isotopes, includi ng C oulomb corrections to 
the nuclear binding energy Il5l.|l6l , but neglecting the ef- 
fects of degenerate nucleons 'Yfl. The combination of this 
NSE with electron capture rates for approximately 200 
nuclei with A = 45-112, which we have determined here 
and in Ref. 2|, was used to compute the rate of elec- 
tron capture on nuclei and the emitted neutrino spec- 
tra as a function of temperature, density and electron 
fraction. This is similar to treatments used in investi- 
gations of electron capture duringstellar evolution Q 
and in thermonuclear supernovae |18| . The rates for the 
inverse neutrino-absorption process are determined from 
the electron capture rates by detailed balance. Due to its 
much smaller |Q|-value, the electron capture rate on the 
free protons is larger than the rates of abundant nuclei 
during the core collapse (fig.nj. However, this is mislead- 
ing as the low entropy keeps the protons significantly less 
abundant than heavy nuclei during the collapse. Fig. |3| 
shows that the reaction rate on nuclei, Rh, dominates the 
one on protons, R p , by roughly an order of magnitude 
throughout the collapse when the composition is consid- 
ered. Only after the bounce shock has formed does Rp 
become higher than Rh, due to the high entropies and 
high temperatures in the shock-heated matter that re- 
sult in a high proton abundance. The obvious conclusion 
is that electron capture on nuclei must be included in 
collapse simulations. 

It is also important to stress that electron capture on 
nuclei and on free protons differ quite noticeably in the 
neutrino spectra they generate. The average neutrino 
energy, (E v ), of the neutrinos emitted by electron cap- 
ture on nuclei, can be obtained dividing the neutrino en- 
ergy loss rate (defined by i Yi£i where Si is the en- 
ergy loss rate by electron capture on nucleus i and Yi 
denotes its number abundance) by the reaction rate for 
electron capture on nuclei, Rh- For the neutrino spec- 
trum we adopt the parametrized form as defined in |19| . 
adjusted to reproduce the average neutrino energy {E v ). 
The neutrino emissivity is then obtained by multiplying 
the NSE-averaged electron capture rate by the neutrino 
spectra. 

Fig. |3 demonstrates that neutrinos from captures on 
nuclei have a mean energy 40-60% less than those pro- 
duced by capture on protons. Although capture on nuclei 
under stellar conditions involves excited states in the par- 
ent and daughter nuclei, it is mainly the larger |Q| -value 
which significantly shifts the energies of the emitted neu- 
trinos to smaller values. Despite that, the total neutrino 
energy loss rate is larger when electron capture on nuclei 
is considered, caused by the increase in the total (nuclei 
plus protons) electron capture rate. The differences in 
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the neutrino spectra strongly influence neutrino-matter 
interactions, which scale with the square of the neutrino 
energy and are essential for collapse simulations. In cur- 
rent simulations 0, |2(| , the low energy portions of the 
neutrino distribution are populated via neutrino-electron 
inelastic scattering of high energy neutrinos produced by 
electron capture on free protons. Electron capture on nu- 
clei produces neutrinos with significantly lower energies, 
accelerating this redistribution process. In this context, 
inelastic neutrino-nucleus scattering, which is usually ig- 
nored, could also be an important process [2ll . 
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FIG. 4: The electron fraction and velocity as functions of the 
enclosed mass at the moment when the center reaches nuclear 
matter densities for a 15 Mq model 0. The thin line is a 
simulation using the Bruenn parameterization Q while the 
thick line is for a simulation using the combined LMP [!j and 
SMMC+RPA rate sets. Both models were calculated with 
Newtonian gravity. 



The effects of this more realistic implementation of 
electron capture on heavy nuclei have been evaluated 
in independent self-consistent neutrino radiation hydro- 
dynamics simulations by the Oak Ridge and Garching 
collaborations [i^, 0] • The basis of these models is de- 
scribed in detail in Refs. and [2£j. Both collapse 
simulations yield qualitatively the same results. Here we 
show a key result obtained by the Oak Ridge collabo- 
ration demonstrating that the effects of this improved 
treatment of nuclear electron capture are twofold. In 
regions close to the center of the star, the additional 
electron capture on heavy nuclei results in more elec- 
tron capture in the new models. In regions where nu- 
clei with A < 65 dominate, the LMP rates result in less 
electron capture. The results of these competing effects 
can be seen in the first panel of Figure which shows 
the distribution of Y e throughout the core when the cen- 
tral density reaches 10 14 g cm -3 , making the transition 
to nuclear matter. The combination of increased elec- 



tron capture in the interior with reduced electron cap- 
ture in the outer regions displaces the velocity minimum, 
which marks the eventual location of shock formation, by 
0.1 Mq. The full effects of these changes on the bounce 
and post-bounce evolution in supernova models will be 
discussed in [ilE^. 

Our calculations clearly show that the many neutron- 
rich nuclei which dominate the nuclear composition 
throughout the collapse of a massive star also dominate 
the rate of electron capture. Astrophysics simulations 
have demonstrated that these rates have a strong impact 
on the core collapse trajectory and the properties of the 
core at bounce. The evaluation of the rates has to rely on 
theory as a direct experimental determination of the rates 
for the relevant stellar conditions (i.e. rather high tem- 
peratures) is currently impossible. Nevertheless it is im- 
portant to experimentally explore the configuration mix- 
ing between pf and sdg shell in extremely neutron-rich 
nuclei as such understanding will guide and severely con- 
strain nuclear models. Such guidance is expected from 
future radioactive ion-beam facilities. 
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